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BACKGROUND 


It  is  well  known  that  the  velocity  response  of  a digital 
moving  target  indicator  (MTI)  is  zero  at  multiples  of  that  doppler 
velocity  known  as  the  blind  velocity  (v^)  which  corresponds  to  the 

radar  pulse  repetition  frequency  (f  ).  The  relationship  also  depends 

on  the  transmitter  frequency  ( f ) and  is  given  by 


_ \ f 3 X 108 

Vb  2 p 2f  T 
t 


(1) 


where  \ is  the  radar  wavelength  and  T is  the  pulse  repetition  interval 
(1/f  ).  Although  the  location  of  the  first  velocity  blind  could  be 

increased  by  reducing  T,  the  minimum  value  of  T is  determined  by  the 

maximum  range  (R  ) at  which  it  is  desirable  to  detect  a target  return, 

m 8 
i.e.,  R = cT/2,  where  c is  the  velocity  of  light  (3  X 10  m/ sec) . 
m 

Consequently,  alternative  procedures  for  extending  the  first  velocity 
blind  have  been  proposed.  These  can  be  subdivided  into  block  stagger 
[1,  2]  and  pulse-to-pulse  stagger  [3-6]  techniques. 

A.  Block  Stagger 

Block  stagger  implies  that  the  number  of  pulses  (N^) 

transmitted  during  a particular  beam  position  are  subdivided  into  N, 

th  ® 

blocks  of  N pulses  each.  For  the  j block  of  pulses,  it  follows 
P 

from  Equation  (1)  that  a specified  velocity  blind  v,  (1)  can  be  generated 

D 

by  choosing  either  a unique  transmit  frequency  ft(j)  or  an  interpulse 

spacing  T(j).  The  first  technique  is  referred  to  as  frequency  agility, 
the  latter  as  block  stagger.  The  composite  velocity  response  consisting 
of  the  individual  block  responses  has  the  first  true  blind  velocity 

when 


Vb  = MjVb(j)  j = l’2 Nb  (2) 

holds  for  all  N,  blocks  and  M.  is  a unique  integer  for  each  block, 
b j 

The  velocity  response  for  an  MTI  processor  using  two  blocks  of  three 
pulses  is  shown  in  Figure  1 for  T(l)  * 1.0  msec  and  ft(l)  = 5.5  GHz. 

The  staggered  response  is  achieved  by  a 20%  increase  in  either  the 
transmit  frequency  or  the  interpulse  spacing  for  the  second  block. 
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It  is  evident  that  while  the  original  blind  velocity  has  been  extended 
by  a factor  of  five,  the  response  contains  a considerable  number  of 
intervals  with  poor  response,  referred  to  as  velocity  "dims"  versus 
a true  "blind."  In  fact,  157.  of  the  sampled  response  values  are  less 
than  -10  dB  and  347,  below  -4  dB. 

B.  Pulse-to-Pulse  Stagger 

Pulse-to-pulse  stagger  implies  a change  in  the  interpulse 
spacing  AT(i)  between  each  successive  transmitted  pulse,  typically 
alternating  between  increased  and  decreased  intervals.  Sometimes,  the 
intervals  are  specified  in  terms  of  stagger  ratios  ( R(  i)}  , i.e., 

AT(i)  = T , (3) 


where  T is  the  unstaggered  interpulse  spacing  corresponding  to  the 
original  vfa  in  Equation  (1)  and  R is  the  average  stagger  ratio;  i.e., 

1 s 

R - r I R(k)  , (4) 

k=l 

where  S is  the  number  of  unique  stagger  ratios  used  and  S < - 1. 

Consequently,  the  staggered-PRF  velocity  blind  v^  is  given  by 

*b  = Vb  R>  (5) 

An  example  of  pulse-to-pulse  stagger  is  shown  in  Figure  2 for  a set  of 

10  pulses  with  a stagger  set  (11,16,13,17}  which  generates  v,  = 14.25  v 

o b 

where  the  response  has  been  plotted  assuming  v = 1 m/sec. 


C.  MTI  Signal  Processor 

The  purpose  of  the  MTI  signal  processor  is  to  remove  the 
clutter  spectrum  while  minimizing  attenuation  of  the  doppler  signal. 

A block  diagram  of  the  MTI  processor  is  shown  in F igure  3.  Because  the 
return  signal  is  not  phase-coherent  with  the  IF  mixer  frequency  (i.e., 
an  unknown  phase  error  <t>  exists)  , it  is  necessary  to  employ  two  similar 
channels  (I  & Q)  to  process  the  baseband  spectrum.  The  I-channel 
processor  heterodynes  the  IF  signal  to  baseband  using  cos  (w^t) , where 

ojj  is  the  angular  IF  frequency,  and  attenuates  the  clutter  spectrum 

using  an  N - tap  transversal  filter  with  coefficients  (A(i,j)}  . 
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XI  (i,i) 


Figure  3.  Digital  MTI  signal  processor  for  one  block  of  N returns 


The  following  analysis  describes  the  1-channel;  however,  the 
results  apply  equally  well  to  the  Q-channel  if  the  symbol  "I"  in  XI(i,j) 
and  YI(k,j)  is  replaced  by  "Q"  and  the  cos  ( ) operator  in  Equation  (6) 
is  replaced  by  a sin  ( ) operator.  The  sampled  doppler  return  in  the 
I-channel  is  given  by 

Xl(i,j)  = cos  [a  ffc(j)  vd  T(i,j)  + <Kj)]  i = 1,2, . . . ,N 

(6) 

j = 1,2, . . . ,Nb 

where 


a 

= 4jt/3 

X 10^  (sec/m) 

ft(j) 

- f + 
o 

(j-l)Af  (Hz) 

j-1 

N 

T(i,j) 

= I 

2.  AT(x,y)  + 2,  AT(x,j) 

y=l 

X=1  X=1 

Hi) 

= R - 

0 

v T(N  , j)  (rad)  . 

a P 

The  MTI  filter  acts  on  a set  of  Nm  consecutive  doppler-return  samples 
according  to  the  algorithm 


YI(k,j)  = 2.  A(x,j)  xi(k  + N - x,j)  k=l,2,...,K  . (7) 

x=l 


If  the  filter  is  "fixed  window,"  only  one  output  is  permitted  for  each 

block  of  N = N pulses;  i.e.,  K = 1.  However,  if  the  filter  is 
p m 

"moving  window"  then  N > N , and  K = N -N  f 1 outputs  are  available. 

pm'  pm 

th 

The  I & Q channel  outputs  are  combined  to  form  the  k residue 
for  the  jth  block  as 

Y(k,j)  =/YI2  (k,j)  + YQ2  (k,j)  , (8) 

which  is  often  approximated  by  some  linear  combination,  e.g.,  larger 

th 

magnitude  plus  one-half  the  smaller  one.  The  k residue  is  then  added 

to  N - N other  residues  and  averaged  in  an  integrator  to  form  the 

pm  th 

residue  for  the  j block.  The  average  of  such  block  residues  is 
then  defined  as  the  velocity  response  of  the  MTI  processor  H(v^)  for 
that  particular  doppler  velocity. 
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A. 


D.  Improvement  Factor  Considerations 


The  improvement  factor  (I)  by  definition  is  the  difference 
in  decibels  of  the  signal-to-clutter  ratio  (SCR)  after  passing  through 
the  MTI  signal  processor  of  Figure  3 to  the  SCR  prior  to  MTI  filtering, 
i • e # , 


I = SCRQ  - SCRi  (dB). 


(9) 


Equation  (9)  can  be  rearranged  as  an  expression  of  signal  gain  (S)  minus 
clutter  gain  (C)  , i.e.. 


I = S - C (dB) 


(10) 


where 


i-i  /' 

\T  J 


“b  0 


|H(v) 


dv 


V 

= 2 / 


b 2 

|H(v) | C . (v)  dv 


C . (v)  = 


v/2jT  a 


2 2 
exp  [ - v / 2 a ] 


and  |h(v) | is  the  amplitude  response  of  the  N^-tap  transversal  filter 


while  a is  the  standard  deviation  (in  m/sec)  of  the  Gaussian  clutter 
spectrum  at  the  MTI  input.  Actual  signal  and  clutter  spectral  density 
scale  factors  are  irrelevant  since  the  ratios  of  output  over  input  are 
considered  in  Equation  (10).  Although  the  doppler  signal  basically 
exists  at  one  velocity,  the  definition  of  I is  based  on  the  average 
signal  gain  (S)  over  the  entire  velocity  response  of  the  MTI  filter. 


i.e.,  zero  to  v^.  The  contribution  of  S to  I is  negligible  compared  to 


C in  most  practical  MTI  processors,  S typically  being  restricted  to  the 
range  ±1  dB  as  compared  to  -20  to  -50  dB  range  for  C.  Consequently, 

I is,  practically  speaking,  a measure  of  clutter  suppression  (-C) . 
Barton  [ 7]  has  developed  an  expression  for  the  improvement  factor  as 


a function  of  the  N -tap  filter  using  binomial  coefficients.  For 
m 


N =3,  this  becomes 
m 


I = 3 + 40  log1Q  [vb/2no] 


(dB) 


(ID 


A limitation  on  I due  to  pulse  staggering  is  given  by  Shrader  [ 8] 


I = 20  log 


10 


0.66  vb 


(7  - 1)<J 


(dB), 


(12) 


where  is  the  unstaggered  blind  velocity  and  y is  the  ratio  of  the 

maximum  interpulse  spacing  to  the  minimum  spacing  or,  equivalently, 
the  maximum  stagger  ratio  to  minimum  ratio.  For  wideband  clutter 
(v^/cr  < 1000),  I is  limited  to  50  dB  or  less  by  either  formula.  The 

MTI  improvement  is  also  limited  by  the  number  of  bits  in  the  A/D 
converter,  e.g.,  I < 52.9  dB  for  n = 9 bits.  Consequently,  values  of 
I exceeding  50  dB  are  meaningless  for  MTI  signal  processors  and  effort 
should  be  directed  toward  exchanging  excess  improvement  for  better 
overall  frequency  response  [ 9]  . 


II. 


PROGRAM  PPABLK 


Coding  documentation  for  this  MTI  simulator  is  found  in  the 
appendix.  The  main  program  PPABLK,  which  is  reproduced  in  Table  A-l , 
is  designed  to  read  the  input  specifications  regarding  the  choice  of 
MTI  tap  gains  (A(i,j)} , stagger  ratios  (R(i,j)} , and  unstaggered  inter- 
pulse spacing  {Tj  for  the  jt'1  block  of  N pulses  processed.  The 
program  converts  the  set  of  S unique  stagger  ratios  and  T.  into  a set 

of  N interpulse  spacings  (AT(i,j)}  for  the  jth  block  and  normalizes 
the  tap  gain  set  (A(i,j)}  to  have  0 dB  noise  power  gain.  Alternatively, 

th 

the  interpulse  spacings  for  the  j block  can  be  read  instead  of  the 
stagger  ratios  by  setting  Tj  = 1.  These  parameters  plus  values  of 

initial  range  (Rq) , initial  transmit  frequency  (f^),  and  frequency 

agility  (Af)  are  passed  to  subroutine  SIGPRO,  which  is  reproduced  in 
Table  A-2,  for  calculation  of  the  velocity  response.  Successive  calls 
to  SIGPRO  provide  amplitude  values  which  are  used  to  compute  the 
average  signal  gain  over  a user  specified  passband,  signal-gain 
histogram,  velocity  response  in  decibels  (dB)  , and  the  improvement  factor  for 
various  values  of  Gaussian  clutter  standard  deviation  (a)  . The  program  is 
described  in  more  detail  in  the  following  subsections.  Program  para- 
meters are  indicated  in  parenthesis  to  coordinate  discussion  with  the 
coding  presented  in  the  appendix. 

A.  Input  Data 

A flowchart  of  the  input  segment  of  PPABLK  is  found  in 
Figure  4.  The  first  card  is  simply  a descriptive  label  (AL)  which 
identifies  the  particular  run  and  is  printed  with  the  output  data. 

The  second  card  consists  of  eight  parameters  which  describe  the  pulse 
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train  and  the  signal  processor.  These  parameters  are  described  in 
Table  1 and  are  identified  both  in  terms  of  the  program  (card)  symbol 
and  the  symbol  used  in  the  narrative,  e.g.,  the  total  number  of  pulses 
per  beam  position  is  identified  as  NSAM/Nt.  The  third  card  contains 

eleven  parameters  which  control  the  output  displays  and  are  described 
in  Table  2,  At  this  point  a test  is  made  using  parameter  IPLOT  to 
determine  whether  the  line  printer  (=  0)  or  graphics  terminal  (=  1) 
will  be  used  to  display  the  velocity  response.  If  the  terminal  is 
used,  then  two  additional  cards  are  read  which  contain  four  additional 
labels  (PLX,  PLY,  PLA,  PLB) . These  labels  are  used  as  X-  and  Y-axis 
labels  for  the  velocity  response  and  sensitivity  study,  respectively. 
The  number  of  stagger  ratios  per  block  (N)  is  then  computed  according 
to  a test  of  the  parameter  (IPROC)  which  identifies  the  MTI  processor 
either  as  moving  window  (=  0)  or  fixed  window  (=  1) . Each  of  the 
blocks  of  data  requires  two  additional  cards,  one  containing  the 

unstaggered  interpulse  spacing  (TPRF)  and  either  the  set  of  N stagger 

ratios  (R(I,J)}  or  the  interpulse  spacings  (DT(I,J)}  for  the  j block, 
the  second  the  set  of  Nm  multipliers  (A(I,J)} . Further  details  regard- 
ing the  format  of  the  data  parameters  are  contained  in  the  program 
listing. 


B.  Computations 


The  computation  portion  of  PPABLK  is  described  by  the 
flowchart  of  Figure  5.  The  first  computation  is  to  determine  the  average 
stagger  ratio  (AVGR)  which  is  found  from  the  average  of  the  S unique 

stagger  ratios.  The  interpulse  spacings  (DT(I,J))  in  the  jth  block 
are  then  computed  from  the  appropriate  stagger  ratios  {R(I,J)j  in 
accordance  with  Equation  (3).  However,  if  TPRF(J)  = 1,  then  AVGR  = 1 
and  DT(I,J)  = R(I, J)*10E-6.  Normalized  filter  coefficients  are 
obtained  by  dividing  the  set  (A(I,J)}  by  the  square  root  of  the  sum 
square  of  the  taps;  consequently,  the  normalized  coefficients  sum 
square  to  unity  (i.e.,  the  wideband  noise  power  gain)  is  0 dB.  The 
computed  interpulse  spacings  and  coefficients  are  printed  for  each 
block  of  data.  In  addition,  the  original  blind  velocity  (VBLD)  and 
the  new  blind  with  stagger  (VBLDWS)  are  computed  in  accord  with 
Equations  (1)  and  (5),  respectively.  If  actual  time  delays  are  entered, 
then  VBLD  is  computed  using  the  average  of  the  interpulse  spacings 
(TAVG)  for  T,  and  the  value  of  VBLDWS  is  set  to  zero. 


The  next  major  computation  is  the  determination  of  the  average 
signal  gain  (SIGDB)  over  the  user-specified  passband.  This  gain  is 
determined  by  averaging  samples  of  the  velocity  response  taken  every 
v^  m/sec  computed  by  subroutine  SIGPRO  from  vm  to  v^  as  specified  by 

user  input  data  on  the  third  card.  At  this  point,  the  velocity  response 
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TABLE  1.  DATA  ELEMENTS  FOR  SECOND  CARD 


*' 


No. 

Card/Report  Symbol 

Description 

1 

NSAM/Nt 

Number  pulses  per  beam  position  (<  150) 

2 

NTAPS/N 

m 

Number  multiplications  (taps)  MTI  filter  (£  15) 

3 

NBLK/Nb 

Number  pulse  blocks  per  beam  position  (5  10) 

4 

NSTRA/S 

Number  of  unique  pulse  stagger  ratios  (£N) 

5 

IPROC 

Type  of  MTI  window,  fixed  (=  1)  , moving  (=  0) 

6 

DEL/Af 

Fractional  charge  transmitter  frequency 
between  blocks 

7 

FOGHZ/f  x 
o 

. io"9 

Initial  transmitter  frequency  (GHz) 

8 

ROKM/R  x . 
o 

io"3 

Initial  targer  range  (km) 

TABLE  2. 

DATA  ELEMENTS  FOR  THIRD  CARD 

No. 

Card/Report  Symbol 

Description 

1 

AMIN 

Minimum  response-plot  velocity  (dB) 

2 

AMAX 

Maximum  response-plot  velocity  (dB) 

3 

VO/v 

o 

Initial  velocity  for  response  plot  (m/ sec) 

4 

VINC/av 

Velocity  increment  for  response  plot  (m/  sec) 

5 

VUP/v 

U 

Final  velocity  for  response  plot  (m/ sec) 

6 

VMIN/v 

m 

Minimum  velocity  for  histogram  (m/sec) 

7 

VMAX/vf 

Maximum  velocity  for  histogram  (m/sec) 

8 

VHIS/v, 

n 

Velocity  increment  for  histogram  (m/ sec) 

9 

DELSGV/Aa 

Initial  value  and  increment  for  clutter o (m/sec) 

10 

NUMVL/Ni 

Number  of  velocity  points  to  compute 
improvement  factor 

11 

I PLOT 

Output  display;  printer  (=  0) , terminal  (=  1) 

13 


HISTOGRAM 
#,  %/2  dB 


Figure  5.  Computation  and  output  segments  of  program  PPABLK. 


values  (H(NV))  (£  1001)  are  converted  to  decibels  and  a histogram 
calculation  is  performed  to  determine  in  which  2-dB  interval  ranging 
from  < -10  dB  to  > +6  dB  the  sample  belongs.  The  histogram  information, 
including  the  total,  plus  number  and  percentage  in  each  interval  is 
printed.  The  average  signal  gain  is  used  later  in  the  computation  of 
improvement  factor. 

The  velocity  response  is  computed  from  vq  to  Vy  every  Av  m/ sec  by 

a second  call  to  SIGPRO.  The  response  is  converted  to  decibels,  with 
-100  dB  set  as  a lower  bound,  and  printed  or  X-Y  plotted  depending  on 
the  user  parameter  IPLOT.  If  IPLOT  = 0,  the  printout  is  obtained  by 
calling  subroutine  VELPRT  (see  Table  A-3)  , whereas  if  IPLOT  = 1 the 
plot  is  obtained  via  a call  to  subroutine  WRT  (see  Table  A-4)  which 
puts  the  output  in  a temporary  file  for  later  processing  by  TEK  4015 
graphics  terminal.  Every  tenth  response  point  is  also  tabulated  on 
the  line  printer.  The  program  is  currently  limited  to  computing  1001 
data  points  for  the  velocity  response. 

The  final  computation  is  a determination  of  the  improvement  factor 
(FACTOR)  which  is  a measure  of  signal-to-clutter  power  ratio  before 
and  after  processing.  This  calculation  is  performed  in  decibels  in 
accord  with  Equation  (10)  and  involves  a computation  of  the  clutter  gain 
(which  is  really  a loss)  after  MTI  filtering.  The  clutter  gain  is 
computed  by  a numerical  approximation  to  Equation  (10)  using  N^  velocity 

points  equally  spaced  between  zero  and  5a  m/sec.  The  improvement  factor 
is  computed  for  20  values  of  clutter  standard  deviation  (a)  equally 
spaced  between  A a and  20  Ao . The  resulting  sensitivity- study  data  are 
tabulated,  and  if  IPLOT  = 1 is  also  stored  for  plotting  on  the  TEK 
4015. 


C.  Subroutine  SIGPRO 

The  subroutine  simulates  the  operation  of  the  MTI  signal 
processor  shown  in  Figure  3 and  described  in  Section  I.C.  The  project  of 
filter  coefficients  (A(i, j)} , normalized  for  0 dB  noise  gain,  and  the 
set  of  interpulse  spacings  (DT(i,j)}  are  passed  from  the  main  program 
via  a COMMON  statement  together  with  the  pertinent  descriptive  para- 
meters read  from  input  data  by  the  main  program.  The  initial  doppler 
velocity  (VO),  velocity  increment  (VDEL) , and  number  of  velocity  points 
(NUMVL)  to  be  computed  are  passed  in  the  subroutine  call.  The  computed 
amplitude-squared  response  (H(NV)}  for  the  NUMVL  velocity  points  is 
returned  to  the  main  program  by  the  COMMON  statement.  The  computations 
are  executed  in  double  precision  to  provide  sufficient  accuracy, 
although  (H(NV)}  is  returned  to  the  main  program  as  a single-precision 
array  of  size  1001  for  processing  by  subroutines  VELPRT  and  WRT. 
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The  flow  chart  for  SIGPRO  is  shown  in  Figure  6 and  the  program  listing 
in  Table  A-2.  The  outer  DO  loop  generates  the  NUMVL  velocity  samples 
and  corresponding  amplitude-squared  response  (H) . The  next  pair  of 
DO  loops  generates  the  individual  I & Q channel  input  samples  (XI  and  XQ) 
for  allNBLK  blocks  of  samples.  The  remaining  nest  of  three  DO  loops 
computes  the  individual  filter  outputs  (YI  and  YQ) , the  resultant 

residue  (RES)  for  the  j*"'1  block,  and  the  final  sum  (T)  of  the  NBLK 
residues. 

D.  Subroutine  VELPRT 

Subroutine  VELPRT  is  used  to  display  the  amplitude 
response  {H(v)}  in  dB  as  a function  of  velocity  using  NPTS  values 
beginning  at  VO  m/sec  in  increments  of  VINC  m/sec.  The  ordinate  is 
limited  to  the  range  AMIN  to  AMAX  in  decibels.  All  aforementioned 
parameters  are  included  in  the  subroutine  call.  The  coding  for  this 
subroutine  is  included  in  Table  A-3. 

E.  Subroutine  WRT 

Subroutine  WRT  is  used  to  write  data  generated  by  PPABLK 
on  a data  file,  designated  as  Unit  3,  in  such  a format  that  it  can  be 
used  by  a plot  routine  on  a TEK  4015  graphics  terminal.  The  call  to 
WRT  includes  two  data  arrays  (X  DATA  & Y DATA),  the  number  of  points 
in  the  arrays  (NPTS),  a curve  designator  (TYPLT)  which  is  sec  to  zeio 
for  the  first  curve  on  a given  plot  or  to  1 designating  other  than  the 
first  curve,  a horizontal-axis  label  (HL)  , a vertical-axis  label  (VL) , 
minimum  and  maximum  values  (XMIN,  XMAX,  YMIN,  YMAX)  for  the  X-axis 
and  Y-axis , respectively , and  the  number  of  curves  (NCURVE)  to  be  drawn 
on  a given  plot.  Before  running  PPABLK  with  the  graphics  option 
(IPLOT  =1.)  the  user  must  state  that  Unit  3 should  be  used  as  a data 
file  and  assign  Unit  3,  i.e.,  the  statements  USE  3,  DATA*FILE."  and 
ASG,A  3."  are  inserted  prior  to  the  @XQT  card  in  the  executive 
control  language  cards.  The  coding  for  subroutine  WRT  is  included 
as  Table  A-4  of  the  appendix, 

F.  Program  PLOT 

The  program  PLOT  is  designed  for  use  with  a TEK  4015 
graphics  terminal.  The  program  reads  data  from  a file  using  the  same 
format  employed  by  subroutine  WRT  to  initiate  the  file.  The  program 
then  calls  a set  of  advanced  graphics  subroutines  [10]  to  display  the 
data  in  graphic  form.  The  graphs  are  labeled  with  the  horizontal  and 
vertical  labels  (HL,  VL)  read  into  PPABLK.  The  axes  are  scaled  according 
to  the  values  (XMIN,  XMAX,  YMIN,  YMAX)  read  from  the  data.  If  more  than 
one  curve  is  drawn  on  a given  plot,  each  curve  will  be  drawn  using  a 
unique  identification  symbol.  After  each  plot  is  drawn  the  program  will 
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return  control  to  the  keyboard  to  allow  time  for  examination  or  to  juake 
a copy.  Depressing  the  return  key  will  clear  the  screen  and  initiate 
the  next  plot.  The  program  will  terminate  when  all  the  data  have  been 
read.  The  user  employs  the  same  control  statements  described  in  the 
aforementioned  WRT  narrative.  A listing  of  Program  PLOT  is  found  in 
Table  A-5 . 


III.  TYPICAL  EXAMPLES 


Several  examples  of  how  to  use  PPABLK  are  presented  to  illu- 
strate its  versatility  and  to  verify  the  simulation  by  comparing  to 
results  presented  in  the  literature.  Except  for  Example  5,  all  plot 
calls  will  be  to  the  line  printer  via  subroutine  VELPRT,  thus  eliminating 
the  two  cards  which  supply  plot  labels  PLX,  PLY  and  PLA,  PLB  used  in 
conjunction  with  the  Tektronix  4015  display.  Consequently,  the  input 
for  the  remaining  examples  consists  of  three  cards  followed  by  two  more 
for  each  block  of  data.  The  first  card  is  simply  the  descriptive 
label  (AL)  which  describes  the  pertinent  features  for  that  design.  The 
next  two  cards  contain  the  parameters  previously  listed  in  Tables  1 
and  2,  and  the  data  values  for  each  example  are  listed  in  the  appropriate 
column  of  Table  3.  If  the  data  parameter  is  the  same  as  the  value  used 
in  the  previous  example,  the  corresponding  entry  is  replaced  by  an  asterisk 
(*) . This  helps  to  pinpoint  the  changes  from  one  example  to  another; 
e.g.,  the  second  example  adds  frequency  agility  to  the  previous  block 
stagger  example  and  the  only  non-asterisk  entry  is  DEL  in  Card  2.  The 
fourth  card  contains  the  unstaggered  pulse  repetition  interval  (TPKF) 
and  the  set  of  stagger  ratios  [R}  , which  are  interpreted  as  correspond- 
ing interpulse  spacings  {DTJ  in  microseconds  if  TPRF  = 1.  The  input 
parameters  for  this  card  are  included  in  the  example  narrative.  Finally, 
the  fifth  card  contains  the  multiplier  coefficients  for  the  I & Q 
channel  digital  filters.  Unless  otherwise  stated  in  the  example,  the 
filter  is  assumed  to  be  the  conventional  three-pulse  canceller  (TPC) 
with  weights  £ 1 , -2,  1)  operating  either  fixed  window  or  moving  window 
according  to  the  value  of  IPROC  in  Card  2 being  1 or  0,  respectively. 

A sensitivity  study  comparison  for  the  five  examples  is  tabulated  in 
Table  4 and  summary  velocity  response  data  is  listed  in  Table  5. 


A.  Example  1 —Block  Stagger 


A block  stagger  design  was  illustrated  in  Figure  1. 

The  data  parameters  used  for  the  second  and  third  cards  are  found  in  the 
first  column  of  Table  3.  The  value  of  TPRF  for  the  first  block  was 
1000  and  1200  psec  for  the  second  block.  As  a consequence,  the  original 
velocity  blind  at  27.3  m/sec  is  increased  by  a factor  of  five.  The 
stagger  ratios  are  {1.,  1.}  for  each  block  since  block,  rather  than 
pulse,  stagger  is  used.  The  histogram  parameters  are  chosen  to  cover 
the  region  to  the  first  blind,  i.e.  VMAX  = 136.,  whereas  the  velocity 
response  is  computed  to  300  m/sec  to  clearly  demonstrate  the  repetition 
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TABLE  3.  DATA  VALUES  FOR  CARDS  2 AND  3 


Card 

1 

2 

3 

4 

3 

No.  2 

NSAM 

6 

* 

10 

* 

3 

NTAPS 

3 

* 

* 

5 

3 

NBLK 

2 

* 

1 

* 

* 

NSTRA 

1 

* 

4 

* 

2 

IPROC 

1 

k 

0 

* 

1 

DEL 

0 

0.01 

0 

* 

* 

FOGHZ 

5.5 

k 

* 

k 

1.3 

ROKM 

5 

k 

* 

* 

* 

No.  3 

-40 

k 

* 

k 

* 

10 

: k 

* 

* 

* 

VO 

0 

k 

* 

* 

* 

VINC 

0.5 

* 

0.05 

2 

3 

VUP 

300 

* 

15 

400 

600 

0 

JU 

* 

* 

k 

VMAX 

136 

* 

14.25 

385 

575 

VHIS 

1 

* 

0.25 

5 

* 

DELSGV 

0.05 

* 

0.0018 

0.05 

* 

NUMVL 

201 

* 

* 

* 

k 

I PLOT 

0 

J. 

* 

* 

l 

*Implies  no  change  from  previous  entry  in  same  row. 


of  the  response  at  multiples  of  v^.  It  is  apparent  from  Table  4 that 

block  stagger  provides  excellent  clutter  attenuation;  however,  it  is 
also  apparent  from  Table  5 that  signal  attenuation  exists  over  a 
wide  portion  of  the  velocity  response. 
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TABLE  4.  SENSITIVITY  STUDY  COMPARISONS 


o(m/ sec) 

System,  I (dB) 

Block  (2) 
(TPC) 
(Fig.  1) 

Block  (2)  / 
1%  Agility 
(Fig.  7) 

Pulse  (10) 
(TPC) 
(Fig.  2) 

Pulse  (10) 
(5-CHEB) 
(Fig.  8) 

Pulse  (3) 
(OPT) 
(Fig.  9) 

0.05 

78.2 

78.1 

56.5 

53.3 

47.0 

0.25 

50.3 

50.1 

42.1 

40.6 

46.2 

0.50 

38.3 

38.1 

35.1 

30.1 

44.4 

1.00 

26.4 

26.3 

26.6 

19.0 

40.5 

TABLE  5.  SIGNAL  HISTOGRAM  STUDIES 


System 

Block  (2) 
(TPC) 
(Fig.  1) 

Block  (2)  / 
17,  Agility 
(Fig.  7) 

Pulse  (10) 
(TPC) 
(Fig.  2) 

Pulse  (10) 
(5-CHEB) 
(Fig.  8) 

Pulse  (3) 
(OPT) 
(Fig.  9) 

Avg  Gain 
(dB) 

-0.8 

-0.9 

-0.4 

-0.3 

-0.0 

> 0 dB 
(%) 

29 

29 

34 

41 

40 

± 2 dB 
(%) 

38 

37 

54 

69 

45 

< -4  dB 
(%) 

34 

34 

12 

4 

22 

B.  Example  2 —Frequency  Agility 

This  example  illustrates  that  according  to  Equation  (1) 
the  blind  velocity  may  be  altered  by  changing  the  transmit  frequency; 
however,  large  changes  are  impractical  and  small  changes  move  the  blind 
by  large  factors.  The  effect  of  frequency  agility  will  be  illustrated 
by  a 1%  change  (DEL  = 0,01)  between  the  two  blocks  used  in  the  first 
example.  Only  the  parameter  DEL  is  changed  as  reflected  in  column  2 
of  Table  3,  Comparison  of  the  data  in  Tables  4 and  3,  and  the  new 
velocity  response  shown  in  Figure  7 with  that  of  Figure  1,  clearly 
demonstrates  the  effect  of  1%  frequency  agility  is  minimal  although  all 
true  blinds  are  eliminated  through  300  m/sec. 
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C.  Example  3 - Pulse-Pulse  Stagger 

The  block  stagger  technique  examples  suffer  from  the 
existence  of  large  amounts  of  attenuation  at  selected  intervals  which 
are  hereafter  referred  to  as  velocity  "dims"  rather  than  true  "blinds." 
An  example  of  pulse-to-pulse  stagger  using  the  stagger  ratios  (11,  16, 
13,  17}  was  shown  in  Figure  2.  Since  the  original  blind  was  selected 
as  1 m/sec,  i.e.,  the  velocity  response  was  normalized  with  respect 
to  the  unstaggered  blind  velocity,  it  follows  from  Equation  (1)  that 
TPRF  = 27273  usee  for  f = 5.5  GHz.  The  data  elements  are  listed  in  the 

third  column  of  Table  3.  The  velocity  response  is  identical  to  that 
provided  by  Figure  44  of  Reference  11.  It  is  apparent  upon  comparison 
of  Figures  1 and  2 or  columns  1 and  3 in  Table  5 that  pulse  stagger  has 
a relatively  uniform  velocity  response.  In  order  to  make  a fair  com- 
parison of  improvement  factors,  the  clutter  increment  was  scaled  by 
1/27  of  the  value  used  in  the  other  examples  to  offset  the  normalized 
blind  velocity  of  1 m/sec.  The  improvement  factor  is  degraded  signifi- 
cantly for  very  narrowband  clutter  relative  to  block  stagger;  however, 
using  50  dB  as  a useful  upper  bound,  the  performance  is  quite  satis- 
factory. 


D.  Example  4 -Higher-Order  Filter 

This  example  illustrates  the  role  of  a five-tap  digital 
filter  which  was  designed  without  stagger  considerations  on  the  velocity 
response.  The  filter  was  originally  designed  by  Houts  and  Burlage  [9] 
to  improve  the  passband  performance  of  an  MTI  filter  operating  at 
f = 5 kHz  and  was  designed  to  remove  narrowband  ground  clutter  (a  = 5 

Hz).  For  purposes  of  closer  comparison  with  other  examples,  the  inter- 
pulse spacing  was  changed  to  1000  psec  (v^  = 27.3  m/sec) . The  stagger 

configuration  is  identical  to  Example  3.  The  set  of  tap  gains 
(-0.28004,  -0.16340,  +0.88868,  -0.16340,  -0.28004}  is  used  in  Card  5. 
Comparison  of  the  velocity  response  shown  in  Figure  8 with  that  of 
Figure  3 or  corresponding  histogram  data  in  Table  5 indicates  that  the 
5-tap  filter  has  a better  passband  performance  achieved  at  the  expense 
of  a 2-7  dB  reduction  in  improvement  factor  as  indicated  in  Table  4. 

E.  Example  5 — Optimum  Tap  Gains 

One  technique  proposed  [12,13]  for  improving  the  clutter 
rejection  for  pulse-to-pulse  stagger  situations  is  to  optimally  select 
the  filter  tap  gains.  One  such  design,  from  a paper  by  Murakami  and 
Johnson  [13]  serves  as  the  final  example.  This  example  also  illustrates 
the  use  of  time  delays  (1000,  1200  usee)  instead  of  stagger  ratios  by 
letting  TPRF  = 1 to  signify  time  delay  interpretations  of  the  stagger 
ratio  entries.  The  set  of  optimized  tap  gains  is  (1.0,  -1.823,  0.833} 
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and  the  other  pertinent  changes  are  listed  in  column  5 of  Table  3.  The 
graphs  obtained  from  the  TEK  4015  are  shown  in  Figures  9 and  10,  the 
latter  illustrating  the  sensitivity  of  the  improvement  factor  to  the 
clutter  deviations.  The  velocity  response  is  essentially  identical  to 
that  presented  in  Figure  9 of  the  aforementioned  paper.  The  labels 
for  the  X and  Y axes  of  these  figures  are  supplied  by  two  additional 
cards  following  Card  3,  with  up  to  50  characters  allowed  to  describe 
the  abscissa  and  30  characters  for  the  ordinate. 

IV.  CONCLUSIONS 

Program  PPABLK  correctly  simulates  MTI  signal  processors 
using  block  versus  pulse  stagger,  fixed  versus  moving  window  filters, 
and  stagger  ratios  versus  interpulse  spacings  description  of  the  pulse 
train.  Five  typical  examples  were  presented  which  illustrated  these 
features  plus  optimal  tap  gain  selection,  frequency  agility,  and  norma- 
lized velocity  response.  The  examples  were  selected  partly  to  verify 
the  proper  performance  of  the  algorithm  and  partly  to  illustrate  the 
effect  on  velocity  response  and  improvement  factor  of  some  of  the 
aforementioned  tradeoffs.  They  were  not  meant  to  be  the  quantitative 
justification  of  one  technique  versus  another,  but  did  illustrate  the 
versatility  of  the  program  and  its  potential  use  to  make  quantitative 
comparisons  of  alternative  design  techniques  for  a particular,  radar 
problem. 


w l 


24 


TARGET  VELOCITY  (m/wc) 
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Appendix.  COMPUTER  LISTINGS  FOR  PPABLK 

The  computer  listings  for  the  simulation  of  an  MTI  signal  pro- 
cessor using  pulse-to-pulse  or  block  are  contained  in  the  five 
tables  found  in  this  appendix.  These  routines  have  been  described  in 
the  main  text  and  further  amplification  is  found  in  the  comment  state- 
ments contained  within  the  program  listings.  The  main  program  PPABLK 
is  described  in  Table  A-l  and  subroutines  SIGPRO,  VELPRT , and  WRT 
are  described  in  Tables  A-2  through  A-4.  Table  A-5  describes  program 
PLOT  which  is  used  in  conjunction  with  the  Tektronix  4015  graphic 
display  terminal.  The  routines  identified  by  call  statements  in  pro- 
gram PLOT  are  described  in  Reference  10. 


J 


♦ 

) 


i 


l,A:  . ; ; i,  VlL  CO'ill 


ftln itrW 


TABLE  A-l.  MAIN  PROGRAM  PPABLK 


C**  PROGRAM  SIMULATES  PIILSE-TO-PULSE  AND/OR  BLOCK  STAGGtH  r-'  T I SY  >TK  . 
C**  INPUT  COEFFICIENTS  A(I,J)  ARE  RESCALED  TC  YILLD  0 Hh  .OlSfc  oaIm. 
C**  SIGNAL  GAIN  IS  MEASURED  EVERY  VHIS  FROM  VMlN  To  V"AX 
C**  SIGNAL  GAIN  IS  HISTOGPAMf’ED  FROM  .LE.  -10  DH  TO  .Gfc.+e>  LH  EACH  ; 
C**  IMPROVEMENT  FACTOR  COMPUTED  FOR  ? 0 VALUES  OF  SIGMAV. 

C**  SEE  FORMAT  40  FOR  ADDITIONAL  DETAILS. 

C«*  EACH  VELOCITY  PLOT  REQUIRES  3 LEAD  CARDS 
C**  1.  LABEL  FOR  PLOTS  <13Ab> 

c**  2.  eight  parameters  (sis,  3F5.2) 

0*  a.  nsam  - total  number  of  samples  i duly  * nbl*.  .ll.  ism 

C**  8.  NTAPS  - NUMBER  OF  FILTER  taps.  ( .LE. IS) 

C**  C.  NtiLK  - iMUMBtP  OF  BLOCKS  (GROUPS).  (.LE.IO) 

c**  D.  nstra  - number  of  unique  stagger  patios  <or  juuKPULsr 

C**  INTERVALS) 

C**  E.  IPROC  - type  OF  PROCESSOR.  MOVING  WINDOW  C:0)  OR  FIXE!  (zl 

C**  F.  DEL  - FREQUENCY  AGILITY  ( FDEL  = DEL  * FO). 

C*»  G.  F06H2  - TRANSMITTER  FREQUENCY  (GHZ) 

C**  . H.  ROKM  - INITIAL  TARGET  RANGE  (KM) 

C**  3.  ELEVEN  PARAMETERS  (4F5.0,  2(F10.0,  F5.01,  F5.0.  2 IS) 

C**  A.  AMIN  - MINIMUM  VELOCITY  RESPONSF  ( DB ) 

C**  B.  AMAX  - MAXIMUM  VELOCITY  RESPONSE  (DR) 

C**  C.  VO  - INITIAL  VELOCITY  VALUE  (M./S) 

C**  D.  VlNC  - VELOCITY  INCREMENT  FOR  RESPONSE  (M/S) 

C**  E.  VUP  - FINAL  VELOCITY  VALUE  (M/S) 

C**  F.  VMIN  - MINIMUM  VELOCITY  FOR  HISTOGRAM  DATA 

C**  G.  VMAX  - MAXIMUM  VELOCITY  FOR  HISTOGRAM  DATA 

C**  H.  VHIS  - HISTOGRAM  INCREMENT  (DEFAULT  = t 'VS ) 

C**  I.  DLLSGV-  CLUTTER  STD.  DEV.  (M/S)  INCREMENT 

C**  J.  NUMVL  - NUMBER  OF  VELOCITY  POINTS  COMPUTED  (.LF.1001)  FOR 

C**  SENSITIVITY  STUDY.  ACCURACY  LOSS  OP  n.2  DM  For  PCI 

C**  PTS  VS.  1001,  with  factor  OF  4 time  SAVINGS. 

C*#  K.  I PLOT  - (=0)  CALL  VELPrT  FOR  LI Nr  PRINTFR  PLOT  OF  VELOCITY 

C*»  RESPONSE  AND  TABLE  OF  SENSITIVITY  STUDY  DATA 

C**  (=1)  CALL  WRT  FOR  GRAPHICS  TERMINAL. 

C**  (1)  VELOCITY  RESPONSE  VO  TO  VHP  tVFKY  Vine 

C**  (2)  SENSITIVITY  STUDY  iwr.  FACTOR  VS.  SIOMAV 

C**  IF  GRAPHICS  TERMINAL  IS  USED  ( IPLOT  = 1)  Then  Two 

C**  ADDITIONAL  LABEL  CARDS  (50A1»30A1)  APE  ItiSLRTLi) 

C**  HERE  FOR  VELOCITY  RESPONSE  aPLX  ( SO  ) ,PL r ( 30 ) A,.N 

C**  CLUTTER  SENSITIVITY  STUDY  (PLA.PLB) 

C**  EACH  VELOCITY  PLOT  REQUIRES  Z * NBLK  ADDITIONAL  CARDS. 

C**  I.  A DESCRIPTION  OF  THE  PULSE-TO-PULSE  STaGGER  In  (LOCK  wO.) 

C**  A.  UNSTaGGEREO  PRF  INTERVAL  IN  MICROSECONDS  (Fb.11) 

t**  B.'A  SET  OF  N STAGGER  RATIOS,  K(I,J),(NFS.O) 

C**  N = MSAM/NbLK  - 1 IF  IPPOC  r 0 . OR 

C**  N = NTAPS  - 1 IF  IPPOC  = 1. 

C**  NOTE.  The  ACTUAL  TIME  DELAYS  IN  MICROSEC  CAI.  Hi  READ  1 14 

C**  PLACE  OF  STAGGER  RATIOS  IF  TPHF(J)  = 1.  IN  PART  A 

C**  2.  A StT  OF  TAP  GAINS  FOR  THE  FIRST  BLOCK  (17X,  7F< i . b/BFrt . S ) 

C**  3.  ADDITIONAL  PAIRS  OF  CARJS  LIKE  1 AND  ? FOR  NBLK- 1 (.LOCK  , 

C************************************ ********  ********** *******♦**♦•»»* 
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TABLE  A-l.  (Continued) 


DIMENSION  R(14,10) ,TPRF(10> . VELC(lOOl) , PL A (SO) . PLB ( JO ) ,PLX ( SO ) .PLY 
S ( 30 ) . PH ( 1 0 ) , N3 ( 1 0 ) ,CLUT (20) ,SIGMAV(20) .FACTOR (20) . /.L(  13) 

COMMON  F0.R0. NTAPS. NDLY.NBLK, FUEL. H( 1001 ) >A( IS. Id) , oT (IS. 1(1) 

SRT2PI  = 2.50o628275 
b «EAU(b.l0.END=10u)  AL 
10  FORMAT  1 13At> ) 

REAL  <b. 20)  NSAM.NTAPS.nbLK.NSTFA, I PROC . DEL . FOGHZ . RnKH 
20  FORMAT ( Sib , 3F5. 2 ) 

REAL.  (b.  30)  AMIN,AMAX,Vu.  VInC,  VUP,  VUIN.  VMAX,VHlS,OLLBGV.NUMy/L  , I PL  JT 
30  FORMAT (4Fb. 0.2 (F10. O.F5.0) ,F5. 0,215) 

FO  = FOGHZ  ♦ I.0E+9 
RO  = ROKM  * 1000. 

TC(iPLOT.tU.O)  GO  TO  34 
RLAulS*33)  PLX.PLY.PLA.PL'I 
33  FORMAT (b0Al,30Al) 

34  WRITE(b»40)  AL.  FO&MZ.  DEL.  NS AM,  ROKM 

40  FORMAT(lHl,//bOX..  PROCESSOR  DYNAMIC  FPEQUENCY  RESPuNSt  * ///  1SX, 

S*  PROCESSOR  CONSISTS  OF  I AND  0 CHANNEL  MTI , DETECTOR  AND  RECOMpIN 
SAT  ION  CIRCUIT,  AND  POST-DETECTION  INTEGRATOR. •//  lsX,  • BOTH  PULSE 
S TO  PULSE  AND/OR  BLOCK  STAGGERING  ARE  PERMITTED  USING  AVERAGE  PRF 
SAND  STAGGLR  RATIO  PARAMETERS. • //  15X.  • FREQUENCY  ABILITY  CAN  PE 

SSIMULATLD  BY  A FREQUENCY  INCREASE  (FDEL  = DFL  X FO I EACH  BLOCK.' 

S //  15X, • UP  TO  15  TAPS/SLOCK  AND  10  BLOCKS  CAN  BE  SIMULATED.  '// 

S 20X.13A6.//'  TRANS.  FRED,  - »,F3.1.»  GHZ,  FREO. AGILITY  = • 

S ,F4 .2.  PULSES  / DWELL  r '.13.'.  INITIAL  TARGET  RANGE  = • 

S ,F3. 1 . • KM . ' //) 

IF ( IPROC  .EQ.  0)  GO  TO  42 
N = NTAPS  - 1 
WRI  TE ( b«  41 ) 

41  FORMAT  (34H0  SIGNAL  PROCESSOR  IS  FIXED  V.'INDOV  ) 

GO  TO  44 

42  N = NSAM/NBLK  - 1 
WRITE (6,43) 

43  FORMAT ( 35H0  SIGNAL  PROCESSOR  IS  MOVING  WINDOW) 

44  NQLY  = N ♦ I 

C**  REAL)  STAGGER  RATIOS  AND  TaP  GAINS,  COMPUTE  INTERPULSf  SPACINGS 
DO  9 J=1.NrLK 

READ(5,bO)TPRF( J) , ( R ( I , J) , 1 = 1 , f I > 

SJ  FOV.aT(I5FS,0> 


REAL.  ( 5. 60  ) (A(I.J).I  =1.  NTAPS) 
bO  F0RMAT(17X,7F9.5  / pFs.5) 

ITt-ST  = IFIX(TPRF  ( J)  ) 

AVGR  = 0. 

DO  1 1=1, NSTRA 

1 A VGR  = A VGR  ♦ R ( I , J ) 

AVGR  = AVtiR  / FLOAT (NSTRA) 

IF ( I TEST .NE . 1 ) GO  TO  7 
TAVG  = AVGR  * 10. **(-6) 

AVGR  = 1. 

7 CONST  = TPRF(J)  ♦ 10.0**(-b)  / AV6R 
DO  2 1=1, N 

2 DT ( 1 , J ) = R ( I , J ) « CONST 
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TABLE  A-l.  (Continued) 


DTlNDLY  ,J)  = TPRF(J)  * 10.O**<-6> 

IF  (ITEST.EQ.l)  OT  ( MOLT  * J ) = TAVG 
WRITE (6*70) J, (0T( I , J) , I=1,NDLY ) 

70  FORMAT (9HOBLOCK  N0*I2.»  - TIME  DELAYS ( SEC ) '.HflO.h  / 31 X . 7F 1 0 . f, ) 
WRI TE ( 6 * 80 ) (A( I , J) . 1 = 1 *NTAPS) 

80  F0RMAT(13X,«  ORIGINAL  TAP  GAINS  ',3Fo.S  / 34X,  7t-9.S/) 

C**  COMPUTE  NORMALISED  FILTER  COEFFICIENTS  WITH  0 D«  NOIc;t  t»A I N 
AVGS  = 0. 

DO  4 I=1,NTAPS 
4 AVGS  = AVGS  + A(I,J>**2 
SRAVGS  = SORT (AVGS) 

DO  8 1 = 1 * NT  APS 
b A ( 1 * J)  = A(I,J)  / SRAVGS 

WRITE (6*90 ) ( A ( I * J ) *1  = 1* NT  APS ) 

90  FORMAT  ( 13X . * NORMALIZED  TaP  GAINS  • * flF  Q . G / 34X.  7>-9. 6/) 

| VBLL  = 0 • IS/ (FOGHZ  * ( 1 . +FLOaT ( J-l ) *DEL) *DT ( NDLY , J 1 ) 

*’  VBLDWS  = VhLD  * AVGR 

IF(ITEST.EO.l)  VBLUWS  = 0. 

9 WRlTE(6*3t>)  VBLD* VbLDWS 

36  FOHMATUHO.IOX,  * ORIGINAL  BLINO  VELOCITY  = '.Fo.!,*  M/S  AND 
% NEW  BLIND  = • *F6.  i • • M/S.M 
FOEL  = DEL  * FO 

C**  COMPUTE  AVERAGE  SIGNAL  GAIN  FROM  VMIN  TO  VM flX  IN  STFPS  OF  VmIS. 
IF(VHIS  .EQ.  0.)  VHIS  = 1. 

NHIS  = 1 ♦ IFIXUVMAX  - VMIN ) / VHIS) 

GAIN  = 0. 

call  sigpro(Vmin,vhis,nhis> 

t**  HISTOGRAM  OF  SIGNAL  GAIN  FROM  -in  TO  +6  DP  / 2 DH  STtPS.  fJHlS  PTG. 

, DO  21  I = 1,10 

21  NB(I>  = 0 

DO  26  NV  = 1 * NHIS 
GAIN  = GAIN  * H ( NV ) 

IF (H( NV ) .GT.  1.0  E-10 ) GO  TO  22 
H(NV ) = -100. 

GO  TO  23 

22  H{NV ) = 10.  » ALOG10(H(NV) ) 

* 23  BL  = - 10. 

IF (H(NV ) .LT.  BL)  NR(1)  = ua<u  ♦ 1 
DO  24  1=  2,9 
BU  = BL  ♦ 2. 

. IF(H(NV) .GE.BL.AND.H(NV) .LT.rU)  NB ( I ) = NR ( I ) + 1 

24  BL  = BU 

26  IF(H(NV).GE.BU)  NH(IO)  = Nu(10)  ♦ 1 
sigain  = gain  / floaT(nhis> 

SIGDB  = 10.  * aLOGIR (SIGAIn) 

WRITE ( 6 * 2B ) VMIN,  V^AX.  ( NB ( I ) * 1 = 1 * 1 0 ) . NHIS 

25  FOHMATdHO,  2UX  * • HISTOGRAM  DATA  FOR  VMIN  = •*F3.n,*  M/S  TO  V'  AX 

*=  dFS.O**  M/S»//»  AMPLITJUE  INTERVAL  (DB)  .LT.-ll)  -10/-H  — M 

S/-6  -6/-4  -4/-2  -2/-0  +0/+2  +2/>4  + 4/*r,  .GE.+o  (ALL  I 

*N  DB)  •//•  NUMBER  OF  SAMPLES /INTERVAL  • , 10 ( 1 5, 3X ) , t TOTAL  = *.I4) 

DO  29  1=1*10 

29  PBII)  = FLOAT ( NB ( I ) * 100 ) / FLOAT(NHIS) 
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TABLE  A-l.  (Concluded) 


I 


WR ITE ( 6, 35 ) (PB (I), 1=1, 10) 

35  FORMAtl  29H0  PERCENTAGE  IN  EACH  INTERVAL  , 10 ( F5. 1 , 3X ) / ) 

<-**  FIND  VELOCITY  RESPONSE  (DB)  FROM  Vu  TO  VUP  IN  STEPS  OF  VINC  (M/S) 
NPTS  = IFIX(l.+( VUP-VO ) /VINC) 

call  sigpro( vo. vinc.npts) 

OO  32  uv  = 1 , NPTS 

IF  t H(NV) .GT.l.OE-lO  ) GO  TO  31 

H(I4V)  - -100. 

GO  TO  32 

31  H ( NV  ) = 10.*  ALOGlO(Hd.V)  ) 

32  CONTINUE 
VtLCU)  = VO 

DO  37  K r 2. NPTS 

37  VELC(K)  = VELC(K-l)  ♦ VINC 
IF(IPLOT.EO.l)  GO  TO  3 t 

CALL  VELPRT ( AMIN » AMAX.NPTS.  VO, VINC.  H) 

GO  TO  39 

38  CALL  WRT(VELC.H, NPTS, O..PLX, PLY, VO. VUP. AMIN. AMAX.l) 

WRITE ( 6, 4o ) ( VELC (K) ,H(K) ,K=1,NPTS,10) 

46  FOHMATUHO.IOX, 'VELOCITY  RESPONSE  DATA  PA IRS-VELC ( VS ). RESPONSE ( OH 
»)».//(2X,8(F6.2* F6.lt 3X) ) ) 

39  CONTINUE 

C**  COMPUTE  IMPROVEMENT  FACTOR  I FOR  SIGMAV  = DELSGV  TO  PO  * DELSGV 
C**  VARY  SIGMAV  FROM  DELSGV  TO  20  * PELSGV 
DO  12  K : 1,20 
CLUT(K)  = 0. 

FACTOR(K)  = 0. 

SIGMAV (K ) = DELSGV  * FLOAT ( K ) 

VULL  = 5.*  SIGMAV (K)  / FLOAT ( NUMVL-1 J 
CALL  SIGPROIO. .VDELtNUMVL) 

C**  COMPUTE  CLUTTER  GAIN  FROM  0.  TO  5*SIGMAV  (M/S)  USING  UUMVL  PTS. 

VEL  = 0. 

DO  11  NV  = l.NUMVL 
VEL  = VEL  ♦ VOEL 

ARG  = 0.5  * (VEL  / SIGMAV(K)  )**2 
CLUT(K)  r CLUT(K)  ♦ H(NV)  * EXP ( -ARG) 

11  CONTINUE 

CLGAIN  = 2.  ♦ VDEL  ♦ CLUT(K)  / (SRT2PI  « SIGMAV(K)  ) 

CLUTDB  = 10.  * ALOGIO(CLGAIN) 

12  FACTOR (K)  = SIGDB  - CLUTDB 
WRITE (6,45 ) NUMVL, SIGDB 

45  FORMATUHO,  5X , * SENSITIVITY  STUDY  DATA  ( NIJVVL  = > , 15,  ' , AVG.  PAS 
SStlAND  SIGNAL  GAIN  = »,F6.3»*  DB)'/) 

WR I TE  ( t>»  55)  (SIGMAV  (K  ) , K=l,20  ) 

55  FORMAT (20H  CLUTTER  SIGMA  (M/S)  , 20F5.2/) 

WRITE(b»65)  (FACTOR (K ) , K=l,20  ) 

65  F QRMAT ( 20H  IMP.  FACTOR  (DB)  , 20F5.1  ) 

IF ( I PLOT  .LQ.  0)  GO  TO  5 

CALL  WRT (SIGMAV,F ACTOR, 20,0. , PL A , PLB , 0 . 0 , SIGMA V ( ?0 ) ,10. ,60. , 1 ) 

GO  TO  5 
100  STOP 
END 
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TABLE  A-2. 


SUBROUTINE  SIGPRO 


SUBROUTINE  SIGPRO ( VC , VbEL , wUMVL ) 

C******** ************  ******************************  * **,*.»*,**.*««*.»*.« 
C**  SIMULATES  I and  0 CHANNEL  MTI  SIGNAL  PROCESSOR  USING  TiaO  NT  APS  ** 
C**  TRANSVERSAL  FILTERS  WITH  RESIDUE  = SGPT ( I**2*0**P) . • ,bLY-NT APS* 1 *♦ 

C**  SUCH  RESIDUES  ARE  SUMMED  FOR  LA  CM  OF  NBLK  BLOCKS  OF  vt  LOL IU  «* 

C*»  RETURNS.  FOR  FIXED-rfINDG*  PROCESSOR,  NDLY  = NT  APS . VtLOCITY  ♦ . 

C**  RESPONSE  = AVERAGE  OF  NBLK*(NDLY-NTaPS+1 ) RESIDUES  Is  C ONTO TEH 
C**  VS.  VELOCITY  (M/S)  FROM  VO  TO  VO  + (NUMVL-1 ) *VDEL.  ♦ * 

c**************** *********  ************************ • ********************* 

DOUBLE  PRECISION  OMEGA .PHASE. PHI , RES, T , XI , xG . Y I . YO 
DIMENSION  XI(I50),XG(1S0) ,RES(10) ,T(10P1 ) 

COMMON  FO »R0, NT APS. NDLY.NULK, FUEL. H( 1001 ) ,A(lS»10).i)T(lb,l0) 

FPIOVC  = 4.191689E-R 
C**  GENERATE  VELOCITY  RESPONSE  H(UV) 

00  100  NV  = l.NUMVL 

V£L  = VO  ♦ FLOAT (NV-1)  * VDEL 

T I = 0. 

T(NV)  = 0 .DO 

C**  6ENERATE  I AND  Q- CHANNEL  VELOCITY-RETURN  SAMPLES 
DO  10'Jsl.NBLK 

XMIT  = (FO  ♦ FLOAT (J-l)  * FDEL)  * FPIOVC 

OMEGA  = XMIT  ♦ VEL 

R r RO  — VEL  * TI 

PHI  = XMIT  • R 

IB  ; J * NOLY 

IA  = IB  - NDLY  ♦ 1 

DO  10  I = IA.IO 

PHASE  = OMEGA  • TI  * PHI 

XI (I)  s DCOS( PHASE) 

XQ(I)  = DSIN ( PHASE ) 

10  TI  : TI  ♦ DT(I-IA*1.J) 

Cm*  FILTERS  I AND  Q CHANNELS  USING  IN-PHASE  AND  QUADRATURE  FILTERS 
DO  30  J : l.NtiLK 
RES(J)  = 0.00 

C**  PROCESSOR  IS  FIXED  WINDOy.  IF  NT  APS  = NDLY. 

DO  25  K s NT APS, NDLY 

Y I s O.DO 

YQ  = O.DO 

DO  20  I r 1 .NT APS 

L = (J-l)  • NDLY  ♦ K - I ♦ 1 

TI  : VI  t A(I.J)  * XI (L ) 

20  YQ  r YQ  ♦ A(I.J)  * XO(L) 

25  RES(J>  = RES ( J ) ♦ bSORI ( YI**2  ♦ YQ**2  ) 
iO  TCNV)  = T(NV)  ♦ RES ( J) 

100  H(NV)  = ( SNGL  ( T(NV ) ) / FLOAT(  NDLK  * ( U(,L  Y - NTAPS  ♦ 1 )>)**<? 
RETURN 
ENO 
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TABLE  A-3.  SUBROUTINE  VELPRT 


SUBROUTINE  VELPRT (Af'IN. AMaX.NUMVL. VO,  VDFL.  M) 
C********************************************************«******»»*»**** 

C**  PRINTS  VELOCITY  RESPONSE  ( H - IN  L)b)  OF  MT I SIGNAL  PROCESSOR  ** 

C**  FROM  VO  TO  VO  + (NUMVL-1 ) *VUEL  (M/S)  IN  INCREMENTS  OF  VDEL . Tut  ** 

C**  RESPONSE  IS  LIMITED  TO  RANGE  AMIN  TO  AMAX  IN  OU.  ** 

C****************  ************************************  *t  ************  t*  * ** 

DIMENSION  H(NUMVL) 

REAL  LINE ( 101 ) 

DATA  DASH»AI.STAR,AX,BLANK/1H-.1hI.1H*,1HX.1H  / 

WRITE  (6,500)  AMIN. AMAX 

500  FORMAT  (lH0.45X.30rlMAGNlTuUE(OH>  VS  VELOCITY  (M/S)  » / / » 7X  » 2Hf)W , OX  . 
$F5. 1 » 95X  »F5. 1 » 4X » 3HM/S ) 

IZERO=IFIX(-AMIN*100./( AMAX-AMIN) )+l 
DO  bio  J=1,NUMVL 
VEL=VO+FLOAT( J-l ) *VDEL 
IF(J.NE.l)  GO  TO  S03 
DO  501  1=1.101 

501  LINE(I)=DASH 
GO  TO  504 

503  DO  502  1=1,101 

502  LINE(I)=BLANK 

504  IFt (IZERO.GT.101) .OP. (1ZERO.LT.1) ) GO  TO  505 
LINE(IZERO)=AI 

505  K=1F IX ( H ( J) *100 ./( AMAX-AMIN ) l+IZERO 
IF(K.GT.lUl)  GO  TO  506 
IF(K.LT.l)  GO  TO  5o7 
LINE(K)=STaR 

GO  TO  510 

506  Krlol 

GO  TO  508 

507  K=1 

508  LINE(K)=AX 

510  WRITE  (6.520)  H ( J ) . LINE , VEL 
520  FORMAT  (1X,F13.3,3X,101A1,F9.2) 

RETURN 

ENU 


TABLE  A-4.  SUBROUTINE  WRT 


SUBROUTINE  WRT ( XDATA  # YU AT  A # NPTS » TYPLT # HL » VL » XMIN»  XMAX  # 
SVMlNr  YMAX  .NCURVE ) 

DIMENSION  XDATa(IO^S) # YUATA ( 10?5) #HL(50) #VL(30> 
IFtTYPLT.EQ.l. ) GO  TO  100 
WRITE ( 3 » 50 ) (HL (I) #1=1# 50) 

WR I TE ( 3 » 70 ) (VL(I) #1  = 1.30) 

50  FORMAT (50A1) 

70  FORMAT ( 3UA1 ) 

WR1  IE  ( 3 # RO  ) NPTS  # NCURVL  # XMIN, X^AX  # Y M I N # Y‘’AX 
80  FORMAT (2Ib#4F10.4) 

WR1 TE ( 3 » 150 ) (XDATA(I).I  = l.NPTS) 

100  WRITE ( 3# 150 ) (YDATA(I).I  = l.NPTS) 

150  FORMAT  (UE13.6) 

RETURN 

END 


t 
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TABLE  A-5.  PROGRAM  PLOT 


DIMENSION  XCATA (1C25 ) » Y D AT  A ( 1C25  J . ML ( 5C ) . VL ( 3 C I 
5C  READ13.  ICC.  END-999)  ( HL  1 1)  *1=1*50  • f VLf  II  • 1=1*30 
ICC  FCRMAT150A1/3CA1) 

READC  3.15C J NPTS.NCURVE .XMIN.XMAX. YMIN.YMAX 
15C  F0RMAT(2I5.4F1C.4) 

F.=  NPT3+1 

F.EAD(3.25C)  (XCAT  A(  I) 

READ ( 3 .25C  ) ( YGATA(  11  .I=2»M> 

XCATA ll)=YC AT A(1)=NPTS 
25C  FORMAT ( 6E13. S ) 

CALL  INTTTC3C) 

CALL  BTNITT 
CALL  TERM(3.1C24) 

CALL  CH?S  I 313) 

CALL  XFRM13) 

CALL  YCRM  l 3 ) 

CALL  XMFRM(l) 

CALL  YMFRMC1) 

CALLV^AME 

CALL  uL IM  X (XM IN  « XKA  X ) 

CALL  CLIMY  ( YMIN'.Yf'AXl 
CALL  TSENG 
CALL  M0VA3F(31C.GCJ 
CC  ICC  I=1.5C 
3CC  CALL  AOUTSTC l.HL( II) 

CALL  MCVA3S( 3Ct GCC I 
CC  3CC  1=1.  30 
CALL  ACUTST ( 1»VL( I ) ) 

CALL  LINEF 
CALL  8AKSP 
3 SC  CONTINUE 

CALL  CHECKtXDATA.YDATA) 

CALL  DCPLAY{XDATA .YDATA ) 

IF ( NCURVC . EG. 1 ) CC  TC  55C 
L=N CURVE-1 
N=NPTS/5 
DC  5CC  K=1»L 

READ( 3 »25C  ) ( YC AT  A ( H . I -2 . M } 

YD AT  A ( 1 ) = NRTS 
CALL  STEFS (N) 

CALL  EYMBLIK) 

CALL  CPLOT(XDATArYDATA) 

5FC  CONTINUE 
5 bC  CALL  TTNPUT(M) 

CC  TC  5C 

999  CALL  FINITTl C.72C) 

STOP 

END 
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